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c3 Magnetic reconnection, a fundamental plasma process associated with a rapid 

dissipation of magnetic energy, is believed to power many disruptive phenomena 
in laboratory plasma devices, the Earth magnetosphere, and the solar corona. 
Traditional reconnection research, geared towards these rather tenuous environ- 
ments, has justifiably ignored the effects of radiation on the reconnection pro- 



ABSTRACT 



re- 
cess. However, in many reconnecting systems in high-energy astrophysics (e.g., 

accretion-disk coronae, relativistic jets, magnetar flares) and, potentially, in pow- 
erful laser plasma and z-pinch experiments, the energy density is so high that 
radiation, in particular radiative cooling, may start to play an important role. 
This observation motivates the development of a theory of high-energy-density 
radiative magnetic reconnection. As a first step towards this goal, we present in 
this paper a simple Sweet-Parker-like theory of non-relativistic resistive-MHD re- 
connection with strong radiative cooling. First, we show how, in the absence of a 
guide magnetic field, intense cooling leads to a strong compression of the plasma 
in the reconnection layer, resulting in a higher reconnection rate. The compres- 
sion ratio and the layer temperature are determined by the balance between ohmic 
heating and radiative cooling. The lower temperature in the radiatively-cooled 
layer leads to a higher Spitzer resistivity and hence to an extra enhancement 
of the reconnection rate. We then apply our general theory to several specific 
astrophysically important radiative processes (bremsstrahlung, cyclotron, and 



inverse-Compton) in the optically thin regime, for both the zero- and strong- 
guide-field cases. We derive specific expressions for key reconnection parameters, 
including the reconnection rate. We also discuss the limitations and conditions 
for applicability of our theory. 

Subject headings: magnetic fields — magnetic reconnection — magnetohydrody- 
namics -- plasmas -- radiation mechanisms: general -- radiation mechanisms: 
thermal 



1. Introduction 

Magnetic reconnection is widely regarded as one of the most important and ubiquitous 
phenomena in plasma physics, with numerous important applications in laboratory plasma 
physics, heliophysics, and astrophysics. Most of the magnetic reconnection research so far has 
been driven by our desire to understand magnetic dissipation in various space, solar, and 
laboratory plasmas, with applications to solar flares, substorms in Earth magnetosphere, 
and sawtooth disruptions in tokamaks. Importantly, all of these environments are relatively 
tenuous, low-energy-density, and optically thin, and are adequately described as a collection 
of non-relativistic charged particles whose numbers are conserved, with no photons. 

Magnetic reconnection has also been frequently invoked in a large variety of astrophys- 
ical contexts, especially in high-energy astrophysics. Examples include accretion disks and 
their coronae and large-scale magnetospheres, jets, gamma-ray bursts (GRBs), pulsar mag- 
netospheres and pulsar winds, flares in soft gamma repeaters (SGRs), etc. Not surprisingly, 
physical insights obtained from solar-, space- and lab reconnection studies have been often 
applied to these astrophysical systems. However, it is important to appreciate that the Uni- 
verse is very diverse and that the range of physical conditions found in various astrophysical 
environments far exceeds that found within our solar system. In particular, there are some 
astrophysical phenomena where, on the one hand reconnection has been hypothesized to 
play an important role and, on the other hand, where the physical parameter regimes are 
really quite different from those in solar flares, Earths magnetosphere, and laboratory plas- 
mas. Therefore, a straightforward extrapolation of the conventional reconnection scalings to 
some of these extreme systems is not justified. The necessity to understand how such sys- 
tems work requires us to develop new theories of magnetic reconnection taking into account 
several physical processes that are not usually included in conventional reconnection studies. 

Some of the most important among these additional physical processes in high-energy 
astrophysical reconnection are those related to the presence of strong magnetic fields and 



hence an overall high level of energy density found in these systems. In particular, at high 
energy densities various radiative processes come into play. In addition, in some magnetically- 
dominated environments with low ambient plasma density, e.g., in pulsar magnetospheres 
and in Active Galactic Nucleus (AGN) jets, one often has to deal with relativistic reconnec- 
tion. And finally, at the most extreme end of high energy density astrophysical reconnection 
- in environments such as magnetar magnetospheres and central engines and inner jets of 
Gamma- Ray Bursts -- the energy density is so high, corresponding to radiation tempera- 
tures in tens of keV, or equivalently, magnetic fields of 10 12 Gauss and more, that, in addition 
to all the radiative processes mentioned above, pair creation should become important inside 
the reconnection layer ( Uzdensky &: MacFadyen|2006 Uzdensky||2008 ) . In the present paper, 
however, we shall concentrate solely on the role of radiative effects in reconnection and will 
leave the relativistic effects and pair creation for a future study. 

In general, radiation may have several important effects on magnetic reconnection dy- 
namics, the most important ones being radiative cooling of the reconnection layer (optically 
thick or optically thin), radiation pressure, and Compton drag (i.e., the radiative resistiv- 
ity due to collisions between electrons and photons, as opposed to electrons and ions as 
in classical Spitzer resistivity). To the best of our knowledge, these aspects of radiative 
magnetic reconnection have not been adequately explored so far, even though they are crit- 
ical for several outstanding problems in modern high-energy astrophysics. There is only a 



handful of references on the subject of reconnection in the presence of radiation (Dorman 



fc Kulsrud|[i"995j |Uzdensky fc MacFadyehl|20~06J |Uzdensky|[2008| |Jaroschek fc Hoshino||2009 



Nalewajko et al. 2010), and clearly much more work needs to be done. Thus, there is a clear 



astrophysical motivation for further effort in this area. 

In addition, we believe that this topic should be of strong interest to experimental High- 
Energy-Density (HED) Physics - - a new exciting branch of modern physics that emerged 
in recent years (e.g., Drake 2006). We therefore anticipate a rapid progress in HED re- 



connection studies facilitated by the advent of new computational and experimental tools 
and capabilities developed in the HED Physics community, including powerful lasers and 
Z pinches. In fact, several HED reconnection experimental studies utilizing laser-produced 



plasmas with mega-gauss magnetic field have already been reported (Nilson et al. 2006 Li 



et al. 2007). 



In this paper we are making the first steps towards building up intuition about the 
role of radiation effects in High-Energy-Density magnetic reconnection. Namely, here we 
concentrate on the simplest case, corresponding to relatively modest energy densities, where 
one needs to include only one of the above-mentioned effects — radiative cooling. In general, 
prompt radiative cooling in both optically thin and optically thick regimes greatly affects the 



energy balance and hence the dynamics of the reconnection layer. Different radiation cooling 
mechanisms may be important in different astrophysical situations, e.g., (1) synchrotron 
and synchrotron-self-Compton in GRB, AGN, and blazar jets; (2) external inverse- Compton 
(EIC) cooling of energetic electrons by powerful ambient soft radiation fields in coronae of 
black holes accreting at a large fraction of the Eddington limit, both in galactic X-ray sources 



and in AGNs (Goodman & Uzdensky 2008); (3) radiation diffusion out of an optically- 
thick pair-dominated reconnection layer in the context of magnetar flares and GRB central 
engines. In general, we believe that any serious effort in this area will require approaching 



the reconnection problem as a radiative-transfer problem (Uzdensky 2008). For simplicity 



however, in this paper we restrict ourselves to the case where the reconnection layer is 
optically thin. We expect this to be a fair assumption at relatively low energy densities, 
(and also for the case of optically thin neutrino cooling in central engines of long GRBs 



and core-collapse supernovae, see, e.g., Kohri et al. 2005). The general optically thick case 



relevant for higher energy densities, is more complicated because of the above-mentioned 
need to solve the radiative transfer problem and is therefore left for a future study. We note 
however, that many of the results derived in this paper, especially in § |2j do not depend on 
the specific properties of the cooling process and therefore are valid for the optically thick 
case as well. 

The outline of our paper is as follows. In § [2] we lay out the general Sweet-Parker-like 
theoretical framework for describing compressible antiparallel magnetic reconnection in the 
presence of strong radiative cooling. In particular, after discussing the main assumptions 
of our model in § 2.1, we work out an analog of the Sweet-Parker model for the case of 



strong plasma compression in § 2.2 with a specific focus on the case of collisional Spitzer 
resistivity (§ |2.3). Then, in § 2.4, we use thermodynamic considerations, namely the heat- 



ing/cooling balance, to elucidate an important connection between radiative cooling and 
the plasma compression inside the layer and obtain general expressions for the reconnection 
layer parameters, including the reconnection rate, in terms of the radiative cooling function. 
In § |3j we apply our general theory to several particularly important radiative mechanisms: 
the bremsstrahlung (§ 3.1) and inverse Compton and cyclotron (§ 3.2). We then devote 



a substantial part of the paper (§ 0| to discussing various conditions of validity of our as- 



sumptions. In particular, in § 4.1 we formulate an evolutionary condition for a reconnecting 
system to be able to reach the strong-compression, strong-cooling regime starting with a 
very hot initial state and we also investigate thermal stability of the reconnection layer. We 
demonstrate that bremsstrahlung does not satisfy the evolutionary and stability conditions 
and thus conclude that, by itself, it cannot lead to the transition to a stable strong-cooling 
regime, but instead may lead to a cooling catastrophe of the reconnection layer, especially 
when aided by another process. In the rest of § [4] we derive the criteria of validity for several 



of our other assumptions and approximations: collisional resistive magnetohydrodynamics 
(MHD) (§ 4.2), the optically thin approximation (§ 4.3), and the Tj = T e assumption (§ 4.4). 



Finally we consider the strong-guide field case in § [5j We summarize our conclusions in § [6j 



2. General Analysis 
2.1. Assumptions 



In this paper we construct a simple Sweet-Parker-like (Sweet 1958 Parker 1957) model 



of magnetic reconnection in the presence of relatively strong radiative cooling. Thus, for sim- 
plicity, we assume a laminar current sheet configuration involving non-relativistic electron- 
ion plasma in a steady state. Strictly speaking, we do not expect these assumptions to be 
satisfied in any real astronomical system with the huge separation of scales (so typical for 
astrophysics) between the global systems size and the thickness of the reconnection current 
layer. We are fully aware that, in reality, even in resistive-MHD without externally driven 
turbulence, very long and thin Sweet-Parker current sheets are likely to be disrupted by the 



secondary tearing instability (Loureiro et al. 2007) and to break up into chains of rapidly 



growing secondary plasmoids ( Samtaney et al.||2009 ; Bhattacharjee et aL]|2009 Cassak et al. 
2009), leading to a very dynamic and non-stationary reconnection process. Correspondingly, 



we do not expect our laminar model to be straight-forwardly applicable to any concrete 
astrophysical situation. However, since our present work represents a pioneering first effort 
laying the foundation for a more realistic radiative reconnection theory in the future, we feel 
that developing a simple Sweet-Parker-like laminar model represents a logical first step. 

Similar to the classical Sweet-Parker analysis, we only consider the two-dimensional 
(2D) problem with no variation of any quantities in the direction of the current. We shall 
mostly focus on the case of purely antiparallel reconnection (with no guide field), leaving the 
more general guide-field case for § [5| In addition, we assume that the plasma is sufficiently 
dense, and hence collisional enough, that resistive MHD with a (smoothly varying) scalar 
resistivity is approximately valid. Thus, we ignore the Hall and other extra terms in the 
generalized Ohm's law, but we do derive the condition for the validity of this assumption. 
We also ignore the effects of viscosity and thermal conduction. 

The key difference between our model and the classical Sweet-Parker model is that we 
do not assume the plasma to be incompressible and instead of the incompressibility condition 
we use the balance between ohmic heating and radiative cooling. In this section we develop 
the general formalism for an arbitrary cooling process and most of the results presented 
here are independent of its exact nature. We specialize our analysis to the case of optically 



thin radiative cooling only at the very end of this section. We then consider several specific 
radiative mechanisms in the next section. 

We would like to remark that a priori (without solving any equations) one should expect 
a large parameter space for the two main assumptions of our model — optically-thin cooling 
and collisional reconnection -- to be valid simultaneously. Indeed, it is generally believed 
that the transition from the collisional to collisionless reconnection takes place when the 
layer thickness becomes comparable to the ion collisionless skin depth, di = c/u P i. If we 
are in a situation where we are approaching this transition, the optical depth across the 
layer is simply r(di) = n e adi, where a is the appropriate scattering cross section. Taking for 
simplicity the Thomson cross-section, o? = (8ir/3) r|, where r e = e 2 /m e c 2 ~ 2.8 x 10 _13 cm 
is the electron classical radius, we get 

r(cfc) ~ n e a T d % = *[™L J = ^ , /™i Kr e 3 ) 1/2 * 100 (n.r 3 ) 1 / 2 , 
6 y m e d e 6 y m e 

where d e = c/u pe is the electron skin-depth. Thus we see that, unless we are dealing with 
extremely huge densities, approaching ~ 10 -4 r~ 3 ~ 5 x 10 33 cm -3 (by which time relativis- 
tic and/or quantum effects would become important anyway), a reconnection layer on the 
brink of the transition to collisionless reconnection is unavoidably optically thin to Thomson 
scattering. 

We shall present and analyze the main equations of our model in the next three sub- 
sections. As mentioned above, we basically follow the classical Sweet-Parker analysis but 
replace the incompressibility condition with the full energy equation. In the spirit of the 
Sweet-Parker model, we can only do rough estimates, valid only up to factors of order unity. 
Nevertheless, we shall sometimes include such factors in our equations, for future reference. 



2.2. Sweet— Parker Reconnection in the strong-compression limit 

We consider a thin reconnection layer (see Fig. [I]) of half-length L and half-thickness 5. 
We denote the direction along the layer as the x direction and the direction across the layer 
as the y direction. We denote the inflow velocity into the layer as t> rec and the outflow 
velocity at the end of the layer as u. We use the subscript "0" to denote the quantities 
measured just upstream of the layer, such as the upstream magnetic field B Q , the upstream 
plasma density n , and the upstream temperature T Q . We also introduce the upstream Alfven 
velocity, 

y/4irn m p 




current 
layer 

Fig. 1. — Reconnection current layer. Magnetic field of strength B (shown by blue lines) is 
reversing across a thin current layer of length 2L and thickness 25 <C 2L. In the presence of 
radiative cooling, most of the plasma thermal energy (resulting from the ohmic dissipation 
of the incoming magnetic energy) is carried away across the layer by photons (shown by 
wiggly red arrows). 



It will also be convenient to define the global upstream Alfven transit time 

L 



tao = 



V A0 



and the corresponding Lundquist number 



«Sn 



LV A0 



>1 



(2) 



(3) 



where rj is the magnetic diffusivity. 

The quantities at the center of the layer will carry no subscripts. We anticipate that 
strong radiative cooling case considered here is characterized by a strong compression of the 



plasma inside the layer (n ^> no), necessary to ensure the pressure balance across it (Dorman 
fc Kulsrud|1995 ). To quantify the degree of this compression, we introduce the dimensionless 
compression ratio 

A=-. (4) 



n 

n 



which, in the regime of interest (strong cooling), will be expected to be a large number, 
A»l. 

The input parameters in our model are the global half-length of the layer, L; the up- 
stream values of the plasma density (no), temperature (To), and the reconnecting magnetic 
field (-Bo); the magnetic diffusivity at the center of the layer, 77; and finally the parameters 
describing the radiative cooling function (to be discussed below). 

Our primary goal is to compute the 5 main output parameters describing the recon- 
nection layer — such as n, T, 5, u, and v rec (or, equivalently the reconnection electric field 
E z = —v rec B /c), in terms of the input parameters. We will therefore need 5 algebraic re- 
lationships between these parameters that we are going to list and discuss in this section. 
In contrast, in the corresponding classical incompressible Sweet-Parker analysis, one would 
have only 4 parameters, because the plasma density inside the layer is fixed to be the same 
as the upstream density by the incompressibility assumption. 

Most, but not all, of these equations will be the same as in the Sweet-Parker model. 
First, we have mass conservation (continuity equation). It is somewhat modified from the 
traditional Sweet-Parker model to take into account the finite compressibility of the plasma. 
Assuming that the plasma density everywhere along the midplane of the layer (including at 
the outflow point N, see Fig. fl| is roughly the same as the density n at the center of the 
layer, we get 

n v Tec L ~ nu5 =>- v rcc L ~ Au5 . (5) 

Next, we consider the set of equations describing the magnetic field: Ohm's law, Fara- 
day's law, and Ampere's law. Here, our treatment is quite standard. From Faraday's law in 
a steady state in two dimensions it follows that the out-of-plane electric field E z is uniform 
across the domain; in the ideal-MHD region just outside the layer (point M), this field can 
be written as E z = —v TCC B /c, whereas at the center of the layer point (O), where both B 
and v vanish, resistive MHD Ohm's law yields E z = r/j z , where r] is the plasma resistivity 
(related to the magnetic diffusivity r] via 77 = 7/c 2 /47r). Estimating j z using Ampere's law as 
j z ~ — cB /4ti5, we get the following important relationship between v rec and 5: 

V ^ V Tec 6 • (6) 

which is the same as in the Sweet-Parker model. 

Furthermore, a quick comparison of the electric field at the inflow point and the ouflow 
point N = (x = L,y = 0), allows us to determine the characteristic reconnected magnetic 
field B\ = B y (L,0) at the end of the reconnection layer: cE z = —v rec B Q = —uB\. (This 
is basically magnetic flux conservation: the amount of the flux that enters the layer per 



unit time equals to the amount of flux that leaves the layer; in other words, magnetic flux 
just reconnects but does not get destroyed). Combining this with the mass conservation 
equation ([5]), we get an estimate 

B 1 = B y (L,0) = P ^ ~ B y A. (7) 

Note that the resulting expression on the right-hand-side differs by a factor of A (which 
we consider to be large) from the corresponding expression in the standard incompressible 
Sweet-Parker theory. As we shall see shortly, this fact will turn out to be significant when 
we determine the outflow velocity at the next step. 

Now, consider the two components of the equation of motion. Since we expect the 
reconnection inflow velocity into the layer, v rec , to be much smaller than the Alfven speed, 
the y component of the equation of motion simply becomes the pressure balance condition. 
Neglecting for simplicity the upstream plasma pressure compared with the magnetic pressure 
(which is justified by for magnetically-dominant coronal environments), we can then write 

P = P(0,0) = S. (8) 

Next, we consider the equation of motion along the layer (at the layer's midplane, y — 0): 

pv x d x v x = - d x P - j z B y /c . (9) 

In the original incompressible Sweet-Parker model, the magnetic tension force (the second 
term on the right-hand side) turns out to be comparable to the pressure-gradient force. 
Then, since one is only interested in getting a rough, order-of-magnitude estimate, one can 
drop the magnetic tension term and integrate the resulting equation from x = to x = L, 
yielding the standard result 

I pu 2 ~ P(0, 0) - P(L, 0) ~ P = ^ . (10) 

where we used the result of the pressure balance across the layer in the last step. This means 
that the outflow velocity in the Sweet-Parker reconnection is the Alfven velocity computed 
with the upstream reconnecting magnetic field Po and the density p (which, we remind the 
reader, is considered to be uniform in the Sweet-Parker model). 

In our situation, however, this is not the case and the tension force turns out to be much 
larger than the pressure gradient force. Indeed, using half the value of reconnecting magnetic 
field B\ at the exhaust point [see eq. (J7|)] to estimate the characteristic reconnected mag- 
netic field inside the layer, we see that the tension force —j z B y /c ~ (B /4:7i5)(B A6/2L) = 
AB 2 /8ttL is by a factor A ^> 1 greater than the pressure gradient force dP/dx ~ B 2 /(87rL). 
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By equating this tension force with the characteristic magnitude of the inertial term, 
pv x d x v x ~ pu 2 /2L, it immediately follows that the outflow velocity at the end of the layer 
is comparable to the upstream Alfven velocity (involving the density and the magnetic field 
just outside of the layer): 

«~ 7= = V M - (11) 

Note that this is different from the usual assumption (which, we claim, is not justified) 
that the outflow velocity in the compressible case is that corresponding to the composite 
Alfven velocity defined with the compressed density inside the layer and the magnetic field 



outside the layer (Parker 1963) - - a result that one would have obtained by taking into 



account the pressure gradient force only. 

Once the compression ratio A is found, solving the rest of the problem (i.e., finding the 
reconnection rate) is rather straight-forward, especially since the remaining equations that 
we need to complete this task [namely, eqs. ([5])-([6]) and (11)] are the same as in the classical 



Sweet-Parker analysis. First, combining Ohm's law (J6J) with the continuity equation (J5l), 



we get v rec L ~ Au5 ~ u (r]/v rec ) A. Then, substituting u ~ Vao from equation (11), we 
immediately arrive at 

Vrec~-^ U J^ A , (12) 

which can be rewritten in terms of the Lundquist number pi) as 

v rcc ~V A vS- 1,2 A 1 ' 2 . (13) 

This means that compression makes reconnection faster by a factor A 1 / 2 , compared with the 
incompressible Sweet-Parker result. Finally, the thickness of the layer is given by 



S „ JL „ L S - x ' 2 A' 1 ' 2 = 6 SP A- 1 ' 2 . (14) 



That is, the compressed layer is thinner than the corresponding incompressible Sweet-Parker 
layer by the same factor of A 1 / 2 . 



2.3. Strong-Compression Reconnection Regime with Spitzer Resistivity 

Next, it is interesting to express the above reconnection layer parameters using the 
actual expression for the resistivity. Assuming that the dominant source of resistivity is 
the perpendicular Spitzer resistivity due to electron-ion collisions (we shall leave the case of 
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Compton-drag resistivity, important in many astrophysical situations, for a future study), 



we can write the corresponding magnetic diffusivity as (Spitzer 1962 Braginskii 1965) 

V± = C v cr e In A 6f'\ 



(15) 



where C v = \/2/3y/ii ~ 0.27 and In A is the Coulomb logarithm. Then, it is straight-forward 
to obtain the following expression for the Lundquist number as a function of density n e and 
temperature 9 e : 

s(n,e e ) = ^ = c-i-±-^e^. (16) 



v± 



What we need in our theory [see eqs. (13)-(14)], is the Lundquist number defined with 
the actual central electron temperature 9 e (which comes in via the resistivity) but with 
the upstream density uq. Furthermore, we would like to express this Lundquist number as 
a function of the compression ratio A. To do this, we make use of the pressure balance 
equation ([8| written as 

P = 2nk B T, (17) 



8tt 



where the factor of 2 reflects the two equal contributions from electrons and ions (whose 



temperatures are assumed here to be the same, see § 4.4). From this, we can express the 



central layer temperature in terms of the compression ratio A: 



k B T{A) 



Bl 



A' 1 k R T, 



B-Lcq ■ 



'11 



lQirn 

where ksT e(l = £>Q/167mo corresponds to the central temperature that the layer would need 
to have without compression (i.e., for A = 1). This temperature can also be written in a 
convenient dimensionless form as 

ksT x Bl \_ mp V|o 

„2 



e e = 



VTLpC'- 



lQirnom e c 2 AA m e c 2 



(19) 



Substituting this expression for the temperature into equation (|16|), we get 
LV A0 



S = S(n , e 



V A0 



%p 



cr 1 — 

v r P In A c 



03/2 



so;^hiA(^) w A °'' 2 - (20) 



Correspondingly, the reconnection velocity (jl 3[) is given by 

^rec 



C 



8C v r e \nA /VaqX- 1 /^e\ 3/4 A s/4 
L \ c J \m n J 



and the the thickness of the layer (|14|) becomes 
L 2 



(21) 



S 2 



C v (Lr e ) In A ±- A-' 6^ 2 = 8C V (Lr e ) In A C^-)~' (^-f 2 A^ . (22) 
Vao V c / \m„/ 
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2.4. Thermodynamics of the Radiativelly-Cooled Reconnection Layer and the 

Compression Ratio 

Finally, we need to consider one more equation that would enable us to determine the 
compression ratio A. To do this, we will need to invoke thermodynamic considerations. As 
we shall see, this part of the problem becomes one-dimensional (ID) in the strong cooling 
limit and effectively decouples from the rest of the analysis. Naturally, one wants to look 
at the energy conservation equation. This equation is usually not needed in the classical 
(incompressible) Sweet-Parker analysis, because there it is replaced by the incompressibility 
condition (although the effects of plasma compressibility have been discussed in the liter- 



ature, see, e.g., Parker 1957, 1963). But since here we are interested in a plasma that is 
strongly compressible, we have to include this equation in our analysis. We note that, even 
though here we are mostly interested in the optically thin case, most of our results in this 
subsection apply, or can be easily extended, to the general optically thick case where one 
has to solve a ID radiative transfer problem. 

In general, energy conservation in a steady state can be viewed as a balance between 
the flux of magnetic enthalpy flowing onto the layer (for simplicity, we here neglect the 
upstream thermal enthalpy compared with the magnetic one; that is, in this paper we restrict 
ourselves to the case of low (3 upstream plasma; however, our analysis can be straight-forwardly 
generalized to the fmite-/3 up stream case) on the one hand and, on the other, the sum of the 
advection of plasma thermal enthalpy and kinetic energy out of the layer, plus (a new 
element!) radiative cooling across the layer. The magnetic enthalpy flux per unit surface area 
is the sum of the inflow of magnetic energy, f rec i?Q/87r, and the work done by the magnetic 
pressure, which is also v rec BQ/8n; thus, the total magnetic enthalpy flux is v rec BQ / 4tt , which of 
course is the same as the Poynting flux, Sp oy nt = cEB /4ii = v tcc Bq/4:Tt (we use the subscript 
"Poynt" to avoid confusion between the Poynting flux and the Lundquist number S). 

Denoting the advective and radiative energy fluxes by -F adv , and -F rad , respectively, we 
can write the energy conservation condition (from the Eulerian viewpoint) as 

Spoynt L ~ — EBq L ~ F adv 5 + F rad L . (23) 

47T 

The total advective energy flux out of the layer can be estimated as 

F adv 5 ~ u {pu 2 /2 + 5nk B T) 5 , (24) 

where the first term is the flux of kinetic energy and the second is the flux of thermal 
enthalpy of a two-species gas. As we showed above, the condition of pressure balance across 
the layer dictates that 2n&sT be equal, or at least comparable, to Bq/8tt. At the same 
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time, the equation of motion along the layer dictates that the plasma kinetic energy is 

pu 2 /2 ~ AB 2 /8tt. Thus, although the two energy fluxes are comparable in the incompressible 

Sweet-Parker case, in the strong-compression A ^> 1 case, the kinetic energy flux always 

dominates over thermal enthalpy flux . It is then easy to see that in this case the ratio of 

the total advective energy flux out of the layer to the total Poynting flux into the layer can 

be estimated as 

F adv o~ ( uA5 



»~>Poynt L \V rec L 

Using the continuity equation ([5]), we thus immediately see that 

-F a dv<5 



(25) 



>->Poynt b 



0(1). (26) 



This is a noteworthy result. It means that, independent of the strength of radiative cooling 
and hence of plasma compression, a finite fraction of the incoming magnetic energy is not 
dissipated into heat (which is then promptly radiated away), but instead is converted directly 
into the mechanical kinetic energy of the outflow. This implies that, strictly speaking, one 
cannot use the total energy conservation to estimate the compression ratio A by simply 
equating the Poynting flux with the radiative losses! What one should do instead, is to 
look at the fate of plasma entropy, i.e., at the balance between the plasma heating due to 
the dissipation of magnetic energy (the ohmic heating) and the plasma cooling both via the 
advection of heat out of the layer and via the radiation losses. In other words, one needs 
to focus on the thermal content of the plasma, as opposed to its total energy. Then one 
can say that the system is in the strong radiative cooling regime if the ohmic heating rate is 
primarily balanced by radiative losses, whereas the heat losses by advection are small. This 
is in contrast with the classical (non-radiative) Sweet-Parker model, where it is straight- 
forward to show that ohmic heating is balanced by the advection of thermal energy out of 
the layer. 

The ohmic heating rate in our model can be estimated as 

,. 2 ,(cB Q \ 2 r\ B 2 v rec Bl 

where rf = Aivq/c 2 is the resistivity, and where we made use of equation (ph. Combining 
this with the mass conservation condition d5J), v Tec L ~ u5A, and with (11), we find 



Qo h m-A-^ = AQ , (28) 

47TT40 

where we defined a characteristic power per unit volume Q = B 2 / \A-ktaq) ■ 
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On the other hand, the heat loss rate due to the advection out of the layer is, roughly, just 
the thermal energy density divided by the advection time along the layer, L/u ~ L/Vao — 

T~A0- 

Qa,dv ~ r~l ~ ~A = Qo ~ Qohm/A , (29) 

L/U 4:71T A0 

where we ignored factors of order unity. 

Thus, not surprisingly, the condition that most of the heat resulting from magnetic 
energy dissipation is promptly radiated away (instead of being advected with the plasma 
outflow) is equivalent to the condition of strong compression (i > 1, u > n ), which is 
what we generally assume in this paper. The exact criterion of when this is the case depends 
on the specific model for radiative cooling; we will consider this in more detail in § |3j 

In general, one can say that the reconnection process partly converts magnetic energy 
into the bulk kinetic energy of the outflow (by performing mechanical work) and partly 
dissipates it into heat (by ohmic heating). In both the Sweet-Parker regime and the radiative 
regime considered here the two parts are roughly comparable. However, whereas in the 
Sweet-Parker case both of these forms of plasma energy are taken out along the layer by 
the plasma outflow, in the strong radiative cooling case the kinetic part of the energy is still 
advected out along the layer but most of the thermal energy is promptly radiated away. 

Provided that we can neglect the advective heat loss compared with the radiative loss, 
we can now write the heating-cooling balance in the optically-thin regime as 

<3rad(™, T) = Qohm ~ A Q , (30) 

where the functional form of the volumetric radiative cooling rate Q ra d on n, and T, (in 
general it may also depend on B, e.g., in the case of cyclotron cooling) depends on the 
specific radiative mechanism. Assuming that this functional dependence is known, this 
equation yields the desired algebraic equation for A: 

Q IBd [n(A),T(A)} ~AQ = A -?°- . (31) 

47rr A0 

Here, the arguments n, T of the radiative cooling function on the left-hand-side are viewed 
as functions of A and of the input plasma parameters n , and B , given by n(A) = n Q A and 



equation (18), respectively. 



It is very important to note that the above equation (31) that determines A does not 
explicitly contain the resistivity 77; it only involves the input plasma parameters no and Bo 
and the length L that are considered to be given. This means that the one-dimensional (ID) 
problem of determining the thermal structure of the layer (i.e., determining the density n 
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and temperature T inside the layer) essentially decouples from the 2D reconnection problem 
itself. This enables us to find the compression ratio A and hence n and T independent of the 
resistivity and the reconnection rate. It is interesting to note that, in general, the resulting 
layer temperature may turn out to be higher or lower than the ambient (upstream) plasma 
temperature. 



3. Specific Radiative Mechanisms 

To make further progress, we need to make a specific choice of the dominant radiative 
cooling mechanism. We shall consider several plausible processes in the following subsections. 



3.1. Bremsstrahlung and Atomic-line Cooling 

In many astrophysical environments, e.g., from the solar corona to the hot gas in galaxy 
clusters, the predominant cooling process is free-free cooling due to binary collisions between 
electrons and ions. Because it is collisional in nature, the cooling rate is proportional to the 
square of the plasma density, i.e., 

Qrad(n,T) = n 2 A(T). (32) 



where A(T) is the so-called cooling function. Correspondingly, equation (31) can be written 
as 

AA(T)^Q n 2 = -^n 2 . (33) 

Atttao 

In general, for a plasma with a non-trivial (e.g., solar) chemical composition, the func- 
tional form of A(T) is quite complicated. In particular, at temperatures between 10 5 K and 
10 7 K cooling is dominated by various strong atomic lines (mostly those of carbon, oxygen, 
and iron) and looks like a mountainous landscape; in solar physics it is usually represented 



by a broken power-law (e.g., Raymond & Smith 1977). However, above about 10 7 K atomic 



lines mostly disappear and the cooling function is dominated by a simple bremsstrahlung 



cooling (Rybicki & Lightman 1979): 



A(T) = kT 1 ' 2 = k Bo A' 1 / 2 = KV AQ {m p /Ak B yi 2 A- l l 2 . (34) 

y/lb7rn k B 

where for a nonrelativistic Maxwellian hydrogen plasma, 

Un 2>i ' ° ^m e k B a fe m ^ 1-4 x 10~ 27 g& erg cm 3 s" 1 K~ 1/2 . (35) 



K ~ 2J1 



3^ 
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where af s = e 2 /he ~ 1/137 is the fine structure constant and g~g is the frequency averaged 
free-free Gaunt factor ( Rybicki fc Lightman|1979 ). Since here we are interested in the regime 
where bremsstrahlung dominates over atomic line emission and where, at the same time, 
the electrons are nonrelativistic, we have to restrict our consideration to the temperature 
range of 10 7 K < T C 5 x 10 9 K; in this range the averaged free-free Gaunt factor varies 



approximately between 1.1 and 1.2 (Karzas & Latter 1961). 



Thus, the bremsstrahlung cooling rate can be expressed as a function of A as 



■\brems / 



Vrad \ A ) 

where 




3^ 



afs r T gs 



T~A0 




3n 



afs r T gs 



lmp ^QoA^ 



8tt 



j r e n L 



(36) 



(37) 



Tt = (TtTIqL 

is the upstream-density Thomson optical depth along the layer. 

Because the above functional form of the cooling function for pure Bremsstrahlung is 



particularly simple, we can solve equation (|33|) corresponding to this case explicitly. We have 

V 2 



A l/2 



'AO 

n L 



K' 1 [k B /m v fl 2 




*T*\ 



2 
AO 



— (afeTr) ~ 
m e 



R 2 
1Q 30 _^_ 



in cgs units . 



(38) 



Then, the condition for the existence of an A ^> 1 solution, corresponding to the strong- 
cooling, strong-compression regime, can be cast as 



V 2 
^»0.5a fs r T 



m P 



m r 



10" 4 T T 



(39) 



or, Bq/uIL > 10 30 in cgs units. 



If this condition is satisfied, we can use equations (21) and (22) to determine the other 
key parameters of the reconnection layer 



^rec 
C 



(3nfHV2C„r- g -^J r -^( V f)' 



m p \ V- 



ttfs T T ) 



-5/2 
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(40) 
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(41) 
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Interestingly, we see that the thickness of the layer becomes independent of the system 
size L in this case; it can be conveniently expressed in terms of the upstream electron 
collisionless skin depth d 2 = m e c 2 / Anrioe 2 = {-Aimorg)' 1 as 



S ,/i llnA c . . 

«e0 V g&Olfs VAe 



where Va £ = B /(4:7inom e ) 1 ^ 2 is the electron Alfven velocity. 



3.2. External Inverse Compton and Cyclotron Cooling 

External Inverse Compton (EIC) and cyclo/synchrotron cooling are two of the most 
important radiative cooling mechanisms in High-Energy Astrophysics, playing a dominant 
role in regulating plasma temperatures in such environments as coronae of accreting black- 
holes and in astrophysical jets, etc. For example, if the plasma is immersed in an isotropic 
external soft radiation field with energy density C/ ra d, then inverse Compton scattering can 
an important cooling mechanism, provided that the characteristic energy of external photons 
is less than the temperature of the layer. This is in fact the case for accretion disk coronae of 
powerful black holes accreting at a sizable fraction of the Eddington limit, e.g., in quasars. 
In these systems, the relatively cool accretion disk provides a powerful source of soft photons 
that very effectively cool the electrons accelerated in coronal reconnection events. Likewise, 
when the plasma has a significant relativistic electron or electron-positron component, as 
is believed to be the case in, e.g., pulsar wind nebulae and in AGN and GRB jets, and is 
immersed in a strong magnetic field, then synchrotron emission may provide the dominant 
cooling mechanism. However, since in the present paper we, for simplicity, focus on non- 
relativistic plasmas where synchrotron emission is absent, we cannot consider this mechanism 
here and we shall therefore leave it to a future study. As for its non-relativistic analog, the 
cyclotron cooling, it turns out that for the non-relativistic zero-guide-field case considered 
in the main part of this paper, it is not effective. The reason for this is the following. By 



virtue of the pressure balance condition (17), the fundamental electron cyclotron frequency, 



Vt ce = eB/m e c , can be expressed in terms of the layer's plasma frequency and the normalized 

1/2 

electron temperature as Vt ce = 2oo pe 6 e ■ Thus, since we assume 9 e -C 1, we see that Vt ce < u pe 
inside the layer; therefore, waves at the electron cyclotron frequency cannot propagate across 
the magnetic field and hence across the layer. There may still be some cooling via higher 
electron cyclotron harmonics, but we neglect it here. We note, however, that in the presence 
of a strong guide field the above pressure balance considerations do not apply and one may 
have a situation where Q ce > u pe and hence where cyclotron cooling can be effective. We 
shall consider this situation in § [5] and devote the rest of this section to inverse- Compton 



cooling. At the same time, we still would like to remark that, even though they are beyond 
the scope of the present paper, the effects of cyclo/syncrotron cooling that one would have 
in the case of a relativistically hot plasma can be treated in a very similar manner. 

Let us investigate the effect of EIC cooling on magnetic reconnection in more detail. 



The optically-thin radiative cooling rate of a single electron can be written as (Rybicki & 



Lightman 1979) 



^;rad '■'■ TT q2 2 



a T cUf5 I 1 \ (43) 

n 3 

where U is the soft radiation energy density, U = U Ta ^. (In the case of cyclotron/synchrotron 

cooling, a similar expression would apply with U being the magnetic energy density, U mag , 

and in the general case when both of these processes are active, U should be the sum of the 

two energy densities, U = C/ ra d + U mag .) 

Next, since in this paper we are focusing on a non-relativistic (,8 < 1) Maxwellian 
plasma with the average electron kinetic energy of < m e c?(3 2 >= 3/2 ksT, we can write the 
radiative cooling rate per unit volume as 

9r>kr,T P 

Q rad = a T cU^- = a T cU— 2 , (44) 

m e c z m e c z 

where we assumed T$ = T e = T to make the last step. Then, using the pressure balance 
condition P = Bq/8ti, we get 

Q rad = °tcU-^- 2 = T -l ^L-^Q . (45) 

Sit m e c z 2 Vao n m e c z 

Interestingly, we see that the volumetric cooling rate is independent of the compression 



ratio A. Then, substituting this expression into equation (31), we can immediately read off 
the compression ratio: 

A ~ ^ rad - Utao - IL JL_ U _ r JiY^L v ^P. 87rU (AQ) 

Q 2m e c 2 2 Vao n m e c 2 4 c m e Bq 



Correspondingly, the normalized layer temperature (19) is 



long: 



n a-1 B #0-1^40 /.«x 

lo7r riom e c z 87r U c 

Then, the condition that A ^> 1 is equivalent to the condition that the layer is sufficiently 

r T »2^^-^ = 4^ — -\ 48 

c U m p Vao o7r(7 
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and the condition that we are dealing with a non-relativistic plasma, 9 e -C 1, as we have 
assumed, becomes 

— J- T T > — . (49) 

By combining these two inequalities, we get the following necessary condition: 

8ttU 




(50) 



Finally, the reconnection layer parameters are: 



Vy 



^rr^^r^i^r, (5i) 



c 2 J V L V c 

*2 



l/2^A0^- 7 / 2 m e /87Tt/y/2 



~ ^^InAr-^j -^J (52) 



4. Additional Validity Conditions 

Let us now discuss the conditions required for the strong-cooling model developed in 
the preceding sections to be valid. 



4.1. The Evolutionary Condition and Thermal Stability of the Reconnection 

Layer 

First, an obvious necessary condition for the strong-cooling, strong-compression regime 



is that equation (31) has a large- A solution, A ~^> 1. However, the actual situation is 
somewhat more subtle. The condition A ^> 1 is just the condition for the existence of a 
stationary strongly-cooled state of the reconnection layer. In addition, however, we must 
impose an extra evolutionary condition for the system to be able to reach this state. As we 
shall see below, this will result in a certain requirement for the radiative cooling function. 
The picture that we have in mind here is the following. The ambient plasma upstream of 
the reconnection layer is rather tenuous; when it just enters the layer, it becomes subject 
to ohmic heating and its temperature rises, whereas its density does not change appreciably 
at first. If radiative cooling can be neglected, one always gets the classical Sweet-Parker 
layer solution, with relatively low density n ~ uq and relatively high temperature T ~ T eq 
(corresponding to A ~ 1). The transition to the strong-cooling, strong-compression A ^> 1 
regime described in the previous sections happens only if that A ~ 1 Sweet-Parker layer 
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becomes unsustainable in the presence of radiative cooling, i.e., if it is able to cool and 
collapse towards the A>1 solution sufficiently rapidly. For this to happen, we must require 
that the radiative cooling of the corresponding A ~ 1 Sweet-Parker solution be stronger 
that the corresponding Ohmic heating, i.e.: 

Qradho, T'eqOo)] > QohmK, T eq (n )] ~ Q , (53) 



where we used equation (28) in the last step. 



Now, assuming that a stationary strong-cooling solution with A ^> 1 does exist, it is 
convenient to make use of the corresponding heating-cooling balance equation Q ra d(A) = 



Qohm(A) ~ AQo [see eq. (30)] and recast the above condition in a form that involves only 



the cooling rate function and the value of A corresponding to the strong-cooling solution: 

QradK, T cq (n )} > Q ohm [n 0l T cq (n )} ~ Q ^ Qohm(A)/A = Q rad (A)/A , (54) 

i.e., 

Q Iad [n(A),T(A)] < AQ vad [n ,T cq (n )} , A > 1 . (55) 

We thus see that the necessary requirement that a reconnecting layer is able to evolve to 
the strong-compression, strong-cooling regime, imposes a certain nontrivial constraint on the 
functional form of the volumetric radiative cooling rate Q(n, T). This means that not every 
radiative process is actually able to cause the system to transition to the strong-cooling 



regime, large- A solution of equation (31), even if this solution exists. For example, it is 



instructive to consider a class of radiative cooling rates having the form 

QT*A(n,T)~ n a T fi . (56) 

This functional form actually reflects many real astrophysically-relevant cooling mechanisms 
such as bremsstrahlung (corresponding to a = 2, (3 = 1/2) and cyclotron/synchrotron and 



inverse-Compton cooling (a = 1, (3 = 1). With this functional form, the condition (55) 
becomes 

a<l + p. (57) 

Thus we see that in the case a = 1, (3 = 1 corresponding to cyclo/synchrotron and ex- 
ternal inverse Compton cooling, important in many high-energy astrophysical environments, 



the condition ( 55 ) is easily satisfied. This suggests that strong radiative cooling of reconnec- 
tion layers by these mechanisms may be important, provided that other necessary conditions 
are also satisfied. 

On the other hand, thermal bremsstrahlung cooling (a = 2, (3 = 1/2) does not satisfy 
the above evolutionary condition. This means that either a large-A solution exists but 
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the system is not able to reach to it because the condition (53) is not satisfied, or the 



condition (53) is satisfied and the reconnection layer starts to cool and compress rapidly but 
a stationary A ^> 1 solution does not exist and so the system experiences a catastrophic 
cooling collapse which can only be stopped by other processes, e.g., by the transition to the 
optically thick regime. (In that latter case, however, it is not clear how one would able to 
avoid a similar cooling catastrophe in the ambient plasma outside the reconnection layer.) 

Moreover, the same conclusion can be generalized to any realistic emission process due 
to binary collisions between electrons and ions, including atomic line cooling, which is the 
predominant cooling process in the solar corona in the temperature range between about 10 5 
and 10 7 K. Because of the collisional nature of these processes, their corresponding cooling 
rate is proportional to the square of the plasma density, Q ra ,d( n , T) = n 2 A(T) , as described 



by equation (32). Hence, a = 2 for any such process, and the condition (57) yields (3 > 1, i.e., 
the cooling function A(T) would have to be rising faster than the temperature. This does not 
appear to be the case for a realistic atomic-line cooling function in the above temperature 
range. 

This conclusion is important, because the fact that the solar corona does not satisfy 
this evolutionary condition is the main reason why prompt radiative cooling is not impor- 
tant in coronal reconnection events (e.g., solar flares). Indeed, one can easily see that a 
solution with large A does exist for solar coronal parameters; therefore, without the above 



evolutionary condition (55), one would expect to see instances of strongly radiatively cooled 
reconnection layers in solar flares, which would have dramatic observational consequences. 
The fact that we don't actually see strong radiative reconnection layers in the solar corona 
is unambiguously explained by the failure of the relevant cooling mechanisms to satisfy the 



condition (55). 



A related very important issue is the thermal stability of the reconnection layer. Let us 
assume that a stationary solution A = A* ^> 1 does indeed exist, Qrad(A.) = Qohm(^4*) = 
A* Q . Then one can easily see that if the function Q r &d(.A)/Q ohra (A) ~ Q ra _ d (A)/A has a 
negative derivative at A = A*, then this solution is stable. If, on the other hand, this function 
has a positive derivative at A — A*, then the solution is unstable: a slight increase in A will 
make Q ra d larger than Qohm, which will lead to a stronger compression so that A will increase 



even further, etc. For the power-law form of the radiative cooling rate, equation (56), we 



have Qrad{A)/Qohm{A) ~ A a " 1 , and hence the condition that the layer is thermally stable 



is a < 1 + f3, which coincides with the above evolutionary condition (57). Thus, EIC leads 



to a stable reconnection layer, whereas the bremssrtahlung solution is unstable. 

It is important to note that in many realistic situations several radiative mechanisms 
are acting simultaneously and, in particular, different cooling processes may dominate at 
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different states of the reconnection layer, which may make the function Q va ^{A)/A non- 
monotonic. For example, it may be possible that one process, characterized by a decreasing 
function Q ra &(A)/A (such as EIC), is responsible for the initial radiative-cooling collapse of 
the hot, low-density, A ~ 1 reconnection layer to a cooler, denser, A ^> 1 state, at which 
point another process, e.g., bremsstrahlung, becomes more effective and takes over as the 
main cooling mechanism. If the cross-over value A c at which Qf^{A^) = Q^ ns (A c ) is 
higher than the value of A corresponding to the stable stationary ElC-cooling solution, see 



equation (46), then that solution is reached before bremsstrahlung becomes effective and the 
above transition to the bremsstrahlung regime does not occur. In the opposite case, how- 
ever, the total (EIC plus bremsstrahlung) cooling rate at the cross-over point A c is higher 
than <5ohm(^4 c ) — A C Q , and the system transitions to the bremsstrahlung cooling regime. 
Then, the catastrophic cooling collapse continues. Eventually, the plasma reaches such low 
temperatures and high densities that either the layer becomes optically thick or recombina- 
tion becomes important and changes the overall physical picture. As follows from the above 
discussion, the condition for this two-stage, hybrid scenario to work can be formulated as 

Q ohm (A c ) ~A C Q < Q™(A C ) = 2 Q^ d c ~ r T -f- — ^ Q , (58) 

Vao n m e c 2 



where we used equation (45) to get the right-hand side, and where A c is the solution of the 



equation Qf^(A) = Q^ 7ls (A). Using expressions (36) and (45), we find that the cross-over 
takes place at 

A = Ac= (^-^ a ^J : ^JL\ V3 . (59) 

\2^2 c \m e n m e c 2 J y ' 

Then, the condition for the transition to the bremsstrahlung-cooling collapse becomes 

U 37T . , 9 , /Vao\ 5 , x 

>— (a is g s )- 2 T T 3 (^) . (60) 



:> " s: ; '-"' ' ' J \ C J 



n m p c 

We believe this two-stage (hybrid) reconnection-layer cooling process may be important 
in some high-energy astrophysical environments, although probably not in the solar corona, 
where the EIC and cyclotron cooling mechanisms are not effective. 



4.2. Collisional Reconnection Condition 

Throughout this paper we have assumed that we are dealing with a resistive-MHD, 
collisional reconnection process. In the absence of a strong guide field, a commonly used 
condition for validity (for non-relativistic electron-ion plasma) of this assumption can be 
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cast as (e.g., Ma & Bhattacharjee 1996 Cassak et al. 2005 Yamada et al. 2006 Uzdensky 
2007a|b|cft 

(61) 



5 S p > di 



where 



di = 
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m e ^4nn r ( 



A-i/2 



(62) 



is the so-called ion collisionless skin depth (r e = e /m e c is the classical electron radius). If 



the condition (61) is not satisfied, then the Sweet-Parker slow resistive reconnection is not 



applicable and one may have a much faster collisionless reconnection regime (attributed to 
either Hall effect or anomolous resistivity). 



Using equation (22), we have 
3 






2 m p Vao 
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(63) 



Thus, the condition that the layer is in the collisional regime, S > di, can be cast as 
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(64) 



More specifically, in the bremsstrahlung case, using equation (38) we see that the re- 
quirement that the layer is collisional reads 
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(65) 



which is easily satisfied for many astrophysical plasmas. 



Similarly, for external inverse Compton cooling, using equation (46), the collisional 
reconnection condition becomes 



l7lU\ 3/2 
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(66) 



4.3. Optically-Thin Condition 

In this paper we have restricted our analysis to the simplest case of an optically thin 
reconnection layer. The reason for this simplification was that we wanted to avoid the 
complexities associated with the full radiative-transfer problem that one has to consider in 
the more general optically thick or intermediate case. 
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The basic necessary condition for the validity of the optically-thin approximation is 
the requirement that the Thomson optical depth across the layer, Tt(S), be small. Using 



equation (22 ), we can express this optical depth in terms of A for the case of Spitzer resistivity 
as 



T T (5)~5a T n~r T A°-~T T A^S-^~T T A^( V -^-) "V^V- /»I^ . (67) 

L V c / \m n / V L 



Because the layer length L is so much larger than r e ~ 3 x lCT 13 cm (by many many 
orders of magnitude) for the absolute majority of potential applications of our theory, we 
see that this optical depth is essentially always much smaller than 1, as we have assumed. 

However, we have to be mindful that Thomson scattering is not the only source of opac- 
ity. In particular, at high densities and low temperatures, free-free absorption, proportional 
to n 2 T~ 7 ^ 2 , may dominate over the electron scattering. In addition, in the relativistic case 
with cyclo/synchrotron cooling, self-absorption may sometimes become important. 



4.4. Electron/Ion Temperature Equilibration Condition 

In this paper we have also implicitly assumed that the electron and ion temperatures are 
equal, or, at least, that the ion temperature is not much higher than the electron temperature. 
If this were not the case, i.e., if we had 7] 3> T e , then our pressure balance condition would 
have been significantly modified; namely, a pressure balance would be established between 
the magnetic pressure outside the layer and the ion pressure inside the layer. As a result, one 
would not have been able to use the pressure balance to determine the electron temperature 
in terms of A, and it is the electron temperature that is needed in order to estimate the 
plasma resistivity and the cooling rate. 

Thus, in order to delineate the realm of validity of our model more accurately, we need 
to figure out when the approximation T e ~ Tj holds. To do this, we follow the approach 



of Uzdensky (2007c). The electron-proton temperature equilibration time teq is by a fac- 
tor m p /m e longer than the electron-ion momentum-transfer collision time T ei ; expressing 
the latter in terms of the electron collisionless skin depth and the collisional magnetic dif- 
fusivity, r ei ~ dl/r], we obtain teq ~ d 2 /rj. It it convenient to recast this estimate using 



expression (14) as follows: 

d 2 d 2 

— ~ Tao A - 

7] 



teq~^ -taqA- 1 ^ (68) 



Because we are interested in the strong-compression regime, /I > 1, we see that the 



25 



condition that electrons and ions are thermally well coupled, teq <C t^ is in fact weaker 
than the condition that the layer is in the collisional regime, 5 > di. 



5. Strong Guide-Field Case 

In the presence of a non-zero guide (i.e., out of plane) magnetic field component B z , the 
situation is somewhat more complicated, in part because one can no longer use the cross- 
layer pressure balance to determine straight-forwardly the central temperature in terms of 
the compression ratio. The guide field provides a certain pressure that resists compression. 
The pressure balance across the layer with an upstream reconnecting magnetic field Bq, 
and an upstream guide field B z q can be written as (again, neglecting the upstream thermal 
pressure) : 

R 2 r2 r2 

8tt 8tt 8tt K J 

where B zc and P c = 2nk B T c are the values of the guide field and the plasma pressure, 
respectively, at the center of the layer. 

Assuming that the guide field compresses in proportion to the gas density, i.e., 

B " " A, (70) 



we can then write the above pressure balance as 

P c = 2An k B T c = ^l + So(l-A 2 ), (71) 

which represents a relationship between T c and A that is somewhat more complicated than 



the simple relationship (18) that we obtained in zero-guide field case. In order to close the 
problem, we need to use the heating/cooling balance to get a second relationship between 
T c and A. In the most interesting for us strong-cooling regime, this balance is basically the 



same as in the zero-guide field case and is given by expression (30): 



Q ra d(An , T c ) = Q ohm ~ AQ . (72) 



Equations (71 ) and (72 ) represent a closed system of two algebraic equations for the un- 
knowns, A and T c , and thus, in principle, completely determine the thermodynamic structure 
of the radiatively-cooled reconnection layer with a finite guide field. 

One particularly simple limit is the case of a strong guide field, B z0 ^> B . In this case, 
the guide field essentially provides a high uniform background pressure and hence effectively 
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enforces incompressibility. Then the pressure balance equation ( 71 ) and the heating/cooling 



balance equation (72) decouple from each other. Namely, in the pressure balance equation 



the gas pressure can be neglected and so the temperature drops out; this equation then only 
involves the compression ratio A and yields a trivial result A = 1 (to lowest order in Bq/ B z q). 



Then, the remaining heating/cooling equation (72) becomes simply 



Qrad(VlO,T c ) = Q 



ohm 



Qc 



Bl 



Attt ao 



(73) 



and one can use it to determine the layer's temperature and therefore the parallel Spitzer 
resistivity tj 



r]±/2 and the corresponding Lundquist number (16): 



S\\ ~2S, 



L V A0 



C v r e In A c 



el' 2 . 



(74) 



Then, the reconnection rate and the thickness of the layer are then given by the usual 
Sweet-Parker scalings: 



V A0 
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L 



s: 
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~2 V A0 L e 



(75) 



For example, in the case of bremsstrahlung, Q rad = n 2 kT 1 ' 2 [see eqs. (32) and (34)- (35) 
we readily obtain 

^tT I fV A o\ 3 m P 



eV 2 




and hence 
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(76) 
(77) 



In the case of EIC cooling, Q rad = IotcV ' n6 e [see eq. (44)], we get 



6 e ~ 



8irU 



Trr 



(78) 



which, interestingly, coincides with the result (47) that we have obtained for the zero-guide 



field case. The main reason for this clearly lies in the fact that, apart from its temperature 



dependence, the EIC radiative cooling rate (44) scales with the compression ratio A in the 
same way (linearly) as the ohmic heating rate Q hm ~ AQ . Therefore, the A factors on 
both sides of the heating/cooling balance cancel out and the resulting temperature turns out 
to be independent of the plasma compression ratio (which is controlled by B z ). 
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In the case of cyclotron cooling, the situation is very similar, we just need to replace U 
in the expression for the radiative cooling rate by the magnetic energy density, U ~ U m3ug 
B^ /8tt. From equation (73) we then immediately get: 



r^ 



n z0 c 



which differs from the corresponding zero-guide-field expression only by a factor of Bq/B^ . 



Recall that, as we discussed in § 3.2 in the B z = case the condition that the electrons 
inside the reconnection layer are non-relativistic (6 e <C 1) coincides with the condition Q ce < 
ujp e (because of the pressure balance and the T e = T t assumption). Then, the cyclotron 
photons cannot propagate effectively through the layer and this is why we chose to neglect the 
cyclotron cooling process for that (relativistic zero-guide field) case. We would like to note, 
however, that in the case of a strong guide field, the restrictions due to the model assumptions 
appear to be not as severe; in particular, there exist a parameter regime where cyclotron 
cooling can be quite effective. This is possible because now the electron temperature is no 
longer determined from the pressure balance and can in fact be quite low. Then, the electron 
cyclotron frequency Q ce can be higher than u pe provided that B^ > Att n m e c 2 . 

We thus see that a more powerful soft radiation field or a stronger magnetic guide field 
lead to a more efficient inverse-Compton or cyclotron cooling, respectively, which results in 
a lower plasma temperature. This, in turn, makes the reconnection process go faster because 
of the T~ 3//2 dependence of the Spitzer resistivity on the temperature. In particular, from 



equation (75), we get 



5 C v r e In A f V AoV ^f 8nU\^ 



'o 



V A0 L V 2 L \ c 
where, once again, U = C/ ra( j + B^ /87i. 



6. Conclusions 

In this paper we presented a theoretical analysis of magnetic reconnection in the pres- 
ence of a strong radiative cooling. Our study was motivated by the increasing interest in 
reconnection in the high-energy astrophysics community and in the laboratory high-energy- 
density community (e.g., laser-plasma and z-pinch communities). Both these areas of re- 
search deal with environments in which the energy density is so high that various radiative 
effects (radiative cooling being just one of them) become important. 
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We would like to stress that this situation can be contrasted with traditional studies 
of magnetic reconnection. These studies are usually motivated by applications to relatively 
tenuous, low-energy-density solar-system environments, such as the solar corona, the Earth 
magnetosphere, laboratory magnetic fusion devices such as tokamaks, and also dedicated 



laboratory reconnection experiments such as the Magnetic Reconnection Experiment (Ya 



mada et al. 1997). One can easily show that in all these systems the effects of radiation 
on the reconnection process are negligible, and this fact justifies the complete omission of 
radiation from the traditional reconnection models. This means, however, that neither any 
concrete results, nor general physical insights from these studies can be applied with any 
degree of confidence to reconnection in high-energy-density environments where radiative 
effects are important. We thus believe that investigating the effects of radiation on magnetic 
reconnection is an important new area of research with potentially significant applications 
to astrophysics. Since this area has been essentially unexplored so far, it should be viewed 
as a new frontier in reconnection research. 

The physics involved in reconnection with radiation is rather complex, as there are 
several different physical effects that radiation can exert, e.g., radiative cooling, radiation 
pressure, radiative resistivity (photon drag), and, in the most extreme astrophysical cases, 
pair creation. We of course do not hope to cover all this complexity in one paper and instead 
take a step-by-step approach. We therefore view the present paper as being a first in a series, 
and so we restricted this study to what we think is the simplest and most fundamental first 
step - - non-relativistic, steady-state, collisional, resistive-MHD reconnection in 2D with 
classical Spitzer resistivity and with radiative cooling. We thus left the above-mentioned 
other potentially important radiative effects for a future study. 

Our analysis in this paper basically followed the classical Sweet-Parker approach but 
added thermodynamic considerations, in particular, the balance between ohmic heating and 
radiative cooling. We devoted most of the paper to the case of strictly anti-parallel reconnec- 
tion, with zero guide magnetic field. We showed that, because of the need to maintain the 
pressure balance with the outside magnetic field pressure, strong radiative cooling inevitably 
leads to a strong compression of the plasma inside the reconnection layer, A = n/n ^> 1, 
where n is the density inside the layer and no is the upstream density just outside the layer. 
We then demonstrated that, in contrast to the incompressible Sweet-Parker case, magnetic 
tension dominates over the pressure gradient as the main force accelerating the plasma out- 
flow from the layer and that the resulting outflow velocity is of order the Alfven speed in the 
uncompressed, upstream plasma. Then, following a Sweet-Parker-like analysis, we derived 
scalings for the the reconnection velocity v rec and the thickness 5 of the layer and showed 
that the former is faster by a factor of A 1 / 2 , and the latter is thinner by the same factor, 
relative to the corresponding classical Sweet-Parker expressions without radiative cooling. 
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We then considered the particular case of Spitzer resistivity and found, not surprisingly, 
that there was another important effect that radiative cooling exerts on the reconnection 
process. Namely, strong cooling makes the layer temperature lower relative to the case 
without cooling, and this, in turn, leads to higher Spitzer resistivity and, correspondingly, a 
lower Lundquist number, which results in faster reconnection. 

We then argued that in order to find the equilibrium layer temperature and the corre- 
sponding compression ratio, one needs to use the condition of balance between the ohmic 
heating and radiative cooling (since in the strong cooling regime the advective heat losses can 
be neglected), and we derived a general equation expressing this balance. The key physical 
idea here is that, whereas in the Sweet-Parker model the matter and the energy go together 
(both advected out along the layer), in the strong radiative cooling case they partly separate: 
while the matter is still advected out of the layer, a significant part of the released energy is 
promptly radiated away across the layer. 

We illustrated these ideas by considering several astrophysically important cooling mech- 
anisms: bremsstrahlung and cyclotron and external inverse Compton emission. For simplic- 
ity, we restricted the analysis in this paper to the optically-thin case, leaving the optically- 
thick and a general intermediate cases to a future study. We derived specific expressions for 
the key reconnection parameters for each of these processes. 

We then considered the restrictions that our various assumptions imposed on the model, 
such as the condition that the layer is in the collisional resistive-MHD regime ( S > di) and 
the condition that it is optically thin. Furthermore, we showed that in addition to the A ^> 1 
strong compression condition, one must impose an evolutionary condition for the system to 
be able to reach the stationary strong-cooling solution starting from its initial hot state. We 
also showed that, related to that condition, there is also a condition of thermal stability of the 
layer. For a cooling rate with a general power-law form, Q ~ ^T 13 , these conditions both 
become simply a < 1 + 0. Interestingly, bremsstrahlung (as well as the solar coronal cooling 
function dominated by various atomic lines) does not satisfy this condition. We believe this 
means that, depending on the particular parameters of a given system, bremsstrahlung is 
either ineffective or leads to a thermal cooling catastrophe. On the other hand, the EIC 
mechanism does satisfy this condition and hence may lead to a stable stationary strong- 
cooling solution. There may also be a hybrid scenario where the initial EIC cooling stage is 
followed by a bremsstrahlung-dominated cooling collapse. 

In the last section of our analysis we discussed the effect of radiative cooling on recon- 
nection with a non-zero guide field, and, in particular, concentrated on the case with a strong 
guide field. We again assumed the Spitzer resistivity and derived specific analytical expres- 
sions for the reconnection rate and the layer thickness. We found that, because a strong 
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guide field effectively enforces incompressibility, there is no acceleration of reconnection due 
to the plasma compression that one has in the zero guide-field case. Thus, the only way in 
which radiative cooling affects (namely, accelerates) the reconnection process in the strong 
guide field case is through its effect on the plasma temperature and then through it -- on 
the Spitzer resistivity and the Lundquist number. 

To sum up, in this paper we made the first strides towards understanding the role of 
radiation (namely, optically thin radiative cooling) in high-energy density magnetic reconnec- 
tion. We hope that, fueled by its potential astrophysical and laboratory high-energy-density 
applications, this important area of research will grow in the coming years. 
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